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Abstract. Let g be a totally positive function of finite type, i.e., g(£) = 
IltLii 1 + ZniSvO" 1 for teK and M > 2. Then the set {e 2 ^ lt g{t - ak) : 
k,l £ Z} is a frame for L 2 (R), if and only if af3 < 1. This result is a first positive 
contribution to a conjecture of I. Daubechies from 1990. So far the complete 
characterization of lattice parameters a, /? that generate a frame has been known 
for only six window functions g. Our main result now provides an uncountable 
class of functions. As a byproduct of the proof method we derive new sampling 
theorems in shift-invariant spaces and obtain the correct Nyquist rate. 



1. Introduction 

The fundamental problem of Gabor analysis is to determine triples (g, a, (3) con- 
sisting of an L 2 -f unction g and lattice parameters a, (3 > 0, such that the set of 
functions Q(g,a,(3) = {e 2n ^ lt g{t - ak) : k, I G Z} constitutes a frame for L 2 (M). 
Thus the fundamental problem is to determine the set (the frame set) 

(1) F{g) = {{a, f3) eR 2 + : Q(g, a, 0) is a frame } . 

It is stunning how little is known about the nature of the set F(g), even after 
twenty years of Gabor analysis. The famous Janssen tie [20] shows that the set 
F{g) can be arbitrarily complicated, even for a "simple" function such as the 
characteristic function g = xi of an interval. 

Under mild conditions, precisely, if g is in the Feichtinger algebra M 1 , then the 
set Fig) is open in [IB] . Furthermore, if g G M 1 , then F{g) contains a neighbor- 
hood U of in Much effort has been spent to improve the analytic estimates 
and make this neighborhood as large as possible [5|[9ll26] . The fundamental density 
theorem asserts that J-{g) is always a subset of {(a, (3) G IR+ : a (3 < 1} [10 ] H H 117] . 
If g G M , then a subtle version of the uncertainty principle, the so-called Balian- 
Low theorem, states that F(g) C {(a, (3) : a/3 < 1} [4,8j. This means that 
{(a, (3) : a(3 < 1} is the maximal set that can occur as a frame set J^ig)- 

Until now, the catalogue of windows g for which F(g) is completely known, 
consists of the following functions: if g is either the Gaussian g(t) = e _7r * 2 , the 
hyperbolic secant g(t) = (e* + e - *)^ 1 , the exponential function e~'*', then J-'(g) = 
{(a, (3) G : a(3 < 1}; if g is the one-sided exponential function g{t) = e~'xR+(£), 
then F{g) = {(a, (3) G : af3 < 1}. In addition, the dilates of these functions 
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and their Fourier transforms, g(t) = (1 + 2ftit)~ l and g(t) = (1 + 4n 2 t 2 )~ l , also 
have the same frame set. The case of the Gaussian was solved independently by 
Lyubarski [25] and Seip [31] in 1990 with methods from complex analysis in re- 
sponse to a conjecture by Daubechies and Grossman [IT] ; the case of the hyperbolic 
secant can be reduced to the Gaussian with a trick of Janssen and Strohmer [22], 
the case of the exponential functions is due to Janssen [19j|2T] . We note that in all 
these cases the necessary density condition a(3 < 1 (or a(3 < 1) is also sufficient 
for Q(g, a, (3) to generate a frame. 

The example of the Gaussian lead Daubechies to conjecture that J-'(g) = {(a, j3) G 
: a(3 < 1} whenever g is a positive function in L 1 with positive Fourier trans- 
form in L l [91 p. 981]. This conjecture was disproved in [18]. 

Surprisingly, no alternatives to Daubechies' conjecture have been formulated 
so far. In this paper we deal with a modification of Daubechies' conjecture and 
prove that the frame set of a large class of functions is indeed the maximal set 
F{g) = {(a,P)eRl:aP<l}. 

This breakthrough is possible by combining ideas from Gabor analysis, approx- 
imation theory and spline theory, and sampling theory. The main observation is 
that all functions above — the Gaussian, the hyperbolic secant, and the exponen- 
tial functions — are totally positive functions. This means that for every two sets 
of increasing real numbers x\ < x 2 < • • • < x^ and yi < y 2 < • • ■ < Vn > N G N, 
the determinant of the matrix [g(xj — Vk)]i<jk<N * s non - ne g a tive. 

Indeed, for a large class of totally positive functions to be defined in (TIB]) we will 
determine the set F{g) completely. Our main result is the following: 

Theorem 1. Assume that g G L 2 (R) is a totally positive function of finite type 
> 2. Then J-{g) = {(a, (3) G : a(3 < 1}. In other words, Q(g,a,/3) is a frame, 
if and only if afi < 1 . 

This theorem increases the number of functions with known frame set from six to 
uncountable. We will see later that the totally positive functions of finite type can 
be parametrized by a countable number of real parameters, see f[T5"j) . Among the 
examples of totally positive functions of finite type are the two-sided exponential 
e~'*' (already known), the truncated power functions g(t) = e~ t t r XR + for r G N, 
the function g(t) = [eT at — e~ bt )xwi+(t) for a, b > 0, or the asymmetric exponential 
g(t) = e at Xm+(—t) + e~ bt Xm+(t), and the convolutions of totally positive functions 
of finite type. In addition the class of g such that J-{g) = {(a, j3) G 1R^_ : a (3 < 
1} is invariant with respect to dilation, time-frequency shifts, and the Fourier 
transform. Since Q(g,a,f3) = Q(g,(3,a), we obtain a complete description of the 
frame set of the Fourier transforms of totally positive functions. For instance, if 
g(t) = (1 + 4vr 2 t 2 )- n for neN, then F(g) = {(a, 0) G R 2 + : a(3 < 1}. 

To compare with Daubechies' original conjecture, we note that every totally pos- 
itive and even function possesses a positive Fourier transform. Theorem [TJ yields 
a large class of functions for which Daubechies' conjecture is indeed true. Fur- 
thermore, Theorem [TJ suggests the modified conjecture that the frame set of every 
continuous totally positive function is J-'(g) = {(a, (3) G M 2 : a(3 < 1}. 
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Our main tool is a generalization of the total positivity to infinite matrices. We 
will show that an infinite matrix of the form [g(xj — Uk)]j P ossesses a left-inverse, 
when g is totally positive and some natural conditions hold for the sequences (xj) 
and (yk) (Theorem [8]). 

The analysis of Gabor frames and the ideas developed in the proof of Theorem [1] 
lead to a surprising progress on another open problem, namely (nonuniform) sam- 
pling in shift-invariant spaces. Fix a generator g G L 2 (M), a step-size h > 0, and 
consider the subspace of L 2 (M) defined by 

V h (g) = {f G L 2 (R) :f = J2 c ^(- ~ kh )} ■ 

For the case h = 1 we write V(g), for short. We assume that the translates g(. — k), 
k G Z, form a Riesz basis for V(g) so that H/H2 x IMh- Shift-invariant spaces 
are used as an attractive substitute of bandlimited functions in signal processing 
to model "almost" bandlimited functions. See the survey [3] for an introduction 
to sampling in shift-invariant spaces. An important problem that is related to the 
analog-digital conversion in signal processing is the derivation of sampling theorems 
for the space V(g). We say that a set of sampling points Xj, ordered linearly as 
Xj < Xj + i, is a set of sampling for Vh(g), if there exist constants A, B > 0, such 
that 

A\\f\\l <£l/(si)| a <S|l/U2 for all / G Vh(g) ■ 

As in the case of Gabor frames there are many qualitative sampling theorems 
for shift-invariant spaces. Typical results require high oversampling rates. They 
state that there exists a 5 > depending on g such that every set with maximum 
gap snpj(xj + i — Xj) = 5 is a set of sampling for V(g) [HEJEI]. In most cases 6 is 
either not specified or too small to be of practical use. The expected result is that 
for V(g) there exists a Nyquist rate and that 6 < 1 is sufficient. And yet, the only 
generators for which the sharp result is known are the 5-splines b n = X[o,i]*' ■ '*X[o,i] 
(n + 1-times). If the maximum gap supJxj+i — Xj) = 5 satisfies 5 < 1, then {xj} 
is a set of sampling for V(b n ) [2]. (This optimal result can also be proved for a 
generalization of splines, the so-called "ripplets" [6]). It has been an open problem 
to identify further classes of shift-invariant spaces for which the optimal sampling 
results hold. 

Here we will prove a similar result for totally positive generators. 

Theorem 2. Let g be a totally positive function of finite type and X = {xj} be a 
set with maximum gap sup^Xj+i — Xj) = 5. If 5 < 1, then X is a set of sampling 
for V(g). 

The theorem will be a corollary of a much more general sampling theorem. 

The paper is organized as follows: In Section 2 we discuss the tool box for the 
Gabor frame problem. We will review some known characterizations of Gabor 
frames and derive some new criteria that are more suitable for our purpose. We 
then recall the main statements about totally positive functions and prove the main 
technical theorem about the existence of a left-inverse of the pre-Gramian matrix. 
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In Section 3 we study Gabor frames and discuss some open problems that are raised 
by our new results. In Section 4 we prove the sampling theorem. 

2. Tools 

2.1. Characterizations of Gabor Frames. There are many results about the 
structure of Gabor frames and numerous characterizations of Gabor frames. In 
principle, one has to check that one of the equivalent conditions for a set Q(g, a, 0) 
to be a frame is satisfied. This task is almost always difficult because it amounts 
to proving the invertibility of an operator or a family of operators on an infinite- 
dimensional Hilbert space. 

In the following we summarize the most important characterizations of Gabor 
frames. These are valid in arbitrary dimension d and for rectangular lattices aZ d x 
(3Z d . We will use the notation MJ = e 2 **/ and T y f = /(• - y), £,t G R d , such 
that 

G{g, a, 13) = {M w T ka g :k,leZ d }. 
Then Q(g, a, 0) is a frame of L 2 (R d ), if there exist constants A,B > 0, such that 

A\\f\\l<J2\(f,M w T ka g)\ 2 <B\\f\\l forall/GL 2 (M d ). 

fc,zez 

If only the right-hand inequality is satisfied, then Q(g,a,(3) is called a Bessel se- 
quence. 

Following the fundamental work of Ron and Shen [26] on shift-invariant systems 
and Gabor frames, we define two families of infinite matrices associated to a given 
window function g G L 2 (R d ) and two lattice parameters a, /3 > 0. 

The pre-Gramian matrix P = P(x) is defined by the entries 

(2) P(x) jk = g(x + ja-%), j, k G Z d . 
The Ron-Shen matrix is G(x) = P(x)*P(x) with the entries 

(3) G(x) kl = '%2g(x + ja-L)g(x + ja-%), k,leZ d . 

j€Z d 

Theorem 3 (Characterizations of Gabor frames). Let g G L 2 {U. d ) and a, /3 > 0. 
Then the following conditions are equivalent: 

(i) The set G(g,a, 0) is a frame for L 2 (R d ) . 

(ii) There exist A, B > such that the spectrum of almost every Ron-Shen matrix 
G(x) is contained in the interval [A,B]: 

o-(G(x)) C [A, B) a.a. xeR d . 

(Hi) There exist A,B > such that 

(4) A||c|| 2 < I °k9(x + ja - |)| 2 < B\\c\\ 2 2 a.a. x G R d , c G £ 2 (Z d ) . 

(iv) There exists a so-called dual window 7, such that ^(7,0,/?) is a Bessel 
sequence and 7 satisfies the biorthogonality condition 

(5) (7, M l/a T k/ pg) = (a/3) d 8 kj0 S li0 , Vh, I G Z d . 
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REMARKS: 1. Condition (iii) is only a simple reformulation of (ii), because (j3J) in 
different notation is just 

\\P(x)c\\ 2 2 = (P( x )c, P(x)c) = (G{x)c,c) x ||c||l. 

2. Condition (iii) as stated has been used successfully by Janssen [19][2T]. The 
construction of a dual window is underlying the proofs that Q(g,a,(3) is a frame 
for the Gaussian g{t) = e _7r * in [TB] and for the one-sided exponential g(t) = 
e-*x R+ (t) [19\. 

3. Condition (iii) establishes a fundamental link between Gabor analysis and 
the theory of sampling in shift-invariant spaces. It says that the set x + aZ is a 
set of sampling for the shift-invariant space Vh(g) with generator g and step-size 
h = 1/(3. Thus Q(g, a, 0) is a frame for L 2 (IR d ), if and only if each set x + aL is a 
set of sampling for V(g) with uniform constants independent of x G M d . 

In our analysis we will use another reformulation of condition (iii). 

Lemma 4. Let g G L 2 (M. d ) and a, (3 > 0. Then the following conditions are 
equivalent: 

(i) The set Q(g,a, 0) is a frame for L 2 {W 1 ) . 

(v) The set of the pre-Gramians {P(x)} is uniformly bounded on i 2 {7L d \ and 
possess a uniformly bounded set of left-inverses, i.e., there exist matrices T(x),x G 
M. d , such that 

(6) T{x)P{x) = I a.a. x G R d , 

(7) \\T(x) \\ op < C a.a. x el 11 . 

In this case, the function 7 defined by j(x + aj) = f3 d T j(x), where x G [0, a) d 
and j G or equivalently 

(8) 7 (ar) = /3 d ^f 0ij (x)x [0 ,^(^-aj), x eR d , 
satisfies the biorthogonality condition (JSJ). 

Proof, (i) =>• (v). If Q(g,a,(3) is a frame, then by Theorem [3](ii) 

^2 I S c k9{x + ja - |)| 2 = (P(x)c,P(x)c) = (G(x)c,c) x ||c|||, 

with bounds independent of x. Consequently G(x) is bounded and invertible on 
£ 2 (Z d ). Therefore the operators P(x) are uniformly bounded on £ 2 (Z d ), and we can 
define T(x) = G(x)~ 1 P*(x). Then 

T(x)P(x) = ((P*(x)P(x))- 1 P*(x))P(x) = Id, 

and 

l|r(x)|| op < iig^-I^iip^iu < a- 1 b 1 / 2 . 

(v) =>• (ii). Conversely, if P(x) possesses a bounded left-inverse T(x), then 
||c|| 2 = ||r(x)P(x)c|| 2 < ||r(x)|| 2 p ||P(x)c|| 2 < C 2 (G(x)c,c) < C7 2 ||P(x)|| 2 p ||c|| 2 , 
and this implies condition (ii) of Theorem [31 
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We next verify that 7 as defined in (JBJ) satisfies the biorthogonality condition (JS} 

(j,M l/a T k/p g) = [ 1 (x)g(x-k/f3)e- 2mlx l a dx 



= f J2l(x + aj)g(x + aj- k/(3) e~ 2nil - x/a dx 
(9) = (3 d [ V T oj (x)g(x + aj - k/(3) e - 27vilx/a dx 

= (3 d [ 5 kfi e- 2mlx/a dx = {a/3)% j0 Si,o • 



The function 7 in (JHJ is a dual window of g, as defined in condition (iv) of 
Theorem [31 provided that G(j, a, j3) is a Bessel sequence. The following result 
gives a sufficient condition. 

Lemma 5. Assume that there exists a (Lebesgue measurable) vector-valued function 
a(x) from IR d to £ 2 (Z d ) with period a, such that 

(10) ^2 °~j( x )9( x + a J ~ I) = 4,o a. a. x eR d . 

If ^2j &Z d sup^g^d |cTj(x)| < oo ; then Q(g,ce,/3) is a frame. Moreover, with 

l( x ) = P d ^2 a i( x )X[o,<x) d ( x - a i)> x e Rd ' 
the set G(j, a, /3) is a dual frame of G(g,a, (3). 

Proof. The computation in shows that 7 satisfies the biorthogonality condition 
flSJ). The additional assumption implies that 

sup \j(x + ak)\ < 00 . 

Consequently, 7 is in the amalgam space W^i 1 ). This property guarantees that 
£(7, a,/3) is a Bessel system [33]. Thus condition (iv) of Theorem [3] is satisfied, 
and Q(g, a, (3) is a frame. The biorthogonality condition (j5J) implies that £(7, a, /3) 
is a dual frame of Q(g, a, (3). 



2.2. Totally Positive Functions. The notion of totally positive functions was 
introduced in 1947 by I. J. Schoenberg [27]. A non-constant measurable function 
g : R — > K is said to be totally positive, if it satisfies the following condition: For 
every two sets of increasing real numbers 



(11) 



X\ < £2 < • • ■ < X N, 



Vi < 2/2 < ■ ■ • < Vn, 



N eN, 
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we have the inequality 

(12) D = det[g(x j -y k )} 1 ^ n >0. 

Schoenberg [28] connected the total positivity to a factorization of the (two-sided) 
Laplace transform of g, 



poo | 

C[g](s) = J j^g(t)dt=: — 



Theorem 6 (Schoenberg [27]). The function g : 1R — > R is totally positive, if and 
only if its two-sided Laplace transform exists in a strip S = {s £ C : a < Res < (3}, 
and 

oo 

with real parameters C, y, 5, 5 U satisfying 

00 

(14) C>0, 7>0, < y + ^Tti* < 00. 



u=l 



A comprehensive study of total positivity is given in the book of Karlin [23]. It 
is known that, if g is totally positive and integrable, then g decays exponentially 
(see [28j p. 340]). In this article, we restrict our attention to the class of totally 
positive functions g G L X (R) with the factorization 

1 M 

as) ^ s)= cm =ce5s ^ l+6us), 

for Mei and real 5 U . This means that the denominator of C[g] has only finitely 
many roots. Equivalently, the Fourier transform of g can be extended to a mero- 
morphic function with a finite number of poles on the imaginary axis and no other 
poles. As noted in [29| p. 247], the exponential factor can be omitted, as it corre- 
sponds to a simple shift of g. In the following we will call a totally positive function 
satisfying ( 1I5|) totally positive of finite type and refer to M as the type of g. 

Schoenberg and Whitney [29] gave a complete characterization of the case when 
the determinant D in ( I12p satisfies D > 0. 

Theorem 7 ( [29]). Let g e L 1 (1R) be a totally-positive function of finite type. 
Furthermore, let m be the number of positive 5 V and n be the number of negative 
b v in f|T5|) . and m + n > 2. For a set of points in ( II ip . the determinant D = 
det[g(xj — yk)]j,k=i,...,N is strictly positive, if and only if 

(16) Xj_ m < yj < x j+n for 1 < j < N. 

Here, we use the convention that Xj = —00, if j < 1, and Xj = oo, if j > N. 

The conditions in f[T6"j) are nowadays called the Schoenberg- Whitney conditions 
for g. They have been used extensively in the analysis of spline interpolation by 
Schoenberg and others (see the monograph [SU]). They will be crucial for our 
construction of a left inverse of the pre-Gramian matrix in (|2]). 
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As a generalization of the pre-Gramians in fl2]), we consider bi- infinite matrices 
of the form 

( 17 ) p =[9(xj-yk)] j>kez , 

where each sequence X = (xj)j e % and Y = (yk)k& C M is strictly increasing. 
Moreover, the sequence (xj)j e z is supposed to be denser in the sense of the following 
condition: 

(a) every interval (yk,yk+i) contains at least one point Xj\ 
.(b) there is an r e N such that \(yk, yk+r) HX\ > r + 1 for all k. 

Our main tool for the study of Gabor frames and sampling theorems will be the 
following technical result. It can be interpreted as a suitable extension of total 
positivity to infinite matrices. 

Theorem 8. Let g G be a totally-positive function of finite type. Letm be the 

number of positive 5 U , n be the number of negative 5 U in ffT5]) . and M = m + n > 1. 
Assume that the sequences (xj)j e z and (y k )kez C M. satisfy condition (C r ). 

Then the matrix P = [g(xj — 2/fe)]- fceZ defines a bounded operator on l^iTD). It 
has an algebraic left-inverse T = [7fcj] fej - 6Z; and 



7fcj = 0, if Xj < y k -rm or Xj > y k+ 



rn ■ 



Proof. We construct a left-inverse T with the desired properties by defining each 
row of r separately. It suffices to consider the row with index k = 0, as the 
construction of all other rows is done in the same way. The goal of the first three 
steps is to select a finite subset of x/s and yk's that satisfy the Schoenberg- Whitney 
conditions. (Our choice of indices is not unique and not symmetric in m and n, it 
minimizes the number of case distinctions.) 

Step 1: Column selection. First, consider the case m,n > and set N := 
(m + n — l)(r + 1), if n > 1, and iV = m(r + 1) + 1, if n = 1. We define an N x N 
submatrix Pq of P in the following way. As columns of Pq, we select columns of P 
between 

fa = -(r + l)m + 1 and k 2 = fa + N - 1 . 

Hence fa = (r + l)(n — 1) for n > 1, and fa = 1 for n = 1. For later purposes, note 
that fa < —m < < n < fa. Therefore, the column with index k = has at least 
m columns to its left and n columns to its right. 

Step 2: Selection of a square matrix. Assumption (C r ) and our definition 
of iV imply that the interval I = (yk 1 + m -i,yk 2 -n+i) contains at least iV points Xj. 
More precisely, for n > 1 we write 

n-2 

{Vkx+m-X, S/fta-n+l) = (V-rm, J/r(n-l)) = [J {Vrv, Vr(u+1)) 

v=—m 

and find at least r + 1 points Xj in each subinterval (y ri/ ,y r r u+ i-\) with — m <v< 
n — 2. This amounts to at least (m + n — 1) (r + 1) = N points in /. If n = 1, we 
have (j/fc 1+fn -i,3/te_n+i) = (y- rm ,yi) and find m(r + 1) points ^ in (y_ rm ,y ) plus 
at least one additional point in (y ,yi). 
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We let 

(19) ji := min{j : xj > y fcl+m _i}, j 2 := max{j : Xj < yk 2 -n+i} ■ 
We have just shown that the set 

X) = {%j ■ ji <j< h) C (y kl+m -i, yk 2 -n+i) 
contains at least X elements. We choose a subset 

that contains precisely X elements and satisfies 

(Vk, Vk+i) n Xq 7^ for fci + m - 1 < k < k 2 - n. 

That is, we choose one point Xj in each interval (y k , 2/fc+i), with ki+m — l<k< 
k 2 — n, and an additional n + m — 1 points Xj G (?/fc 1+m _i, y^-n+i)- Note that 

(20) yhx+m-i < £i = minX < y fel+m < y fe2 _ n < maxXj = £ N < y k2 - n +i- 
Now set 

Vk = Vkn+k-i, 1 <k < N , 

and define the matrix 

p o = - Vk))j,k=i,...,N- 
Then Po is a quadratic X x X-submatrix of P. 

Step 3: Verification of the Schoenberg- Whitney conditions. We next 
show that Po is invertible by checking the Schoenberg- Whitney conditions. First, 
by (EOT) , we have 

& = minXg < y fcl+m = 1m+i' 
By the construction of X' , this inequality progresses from left to right, i.e., 

£j < 2/fei+m-i+j = %+m for \<j<N-m. 
Likewise, we also have 

VN-n = yk 2 -n < maxXg = ^at, 

and this inequality progresses from right to left, i.e., 

7]j < £ j+n for 1 < j < N - n. 

Therefore, the Schoenberg- Whitney conditions (jT6l) are satisfied, and Theorem [7] 
implies that det P > 0. 

Step 4: Linear dependence of the remaining columns of P. We now 

make some important observations. 

Choose indices ko and sGZ with k < k± and m < s < X, and consider the new 
set {rj' k : k — 1, . . . , X} consisting of the points 

yk < yk x < ■ ■ ■ < yk 1+s -2 < yk 1+s <■■■ <y k2 

and the corresponding X x X-matrix Pq = (g(£j — Vk))j,k=i....,N- This matrix is 
obtained from Pq by adding the column (g(£j — yk )) 1< <N as its first column and 



10 



KARLHEINZ GROCHENIG AND JOACHIM STOCKLER 



deleting the column [g(£j — Vs)) 1< <N - Then r] m appears in the m + 1-st column 
of Pq. By (J20J), we see that 

V' m +1 =Vm = Vki+m-1 < 6" 

Consequently, the Schoenberg- Whitney conditions are violated and therefore det Pq = 
by Theorem [3 Since this holds for all s > m, the vector — yk ))i<j<N must 
be in the linear span of the first m columns of Po, namely {g{£,j — Vk))i<j<N for 
k = k%, , . . . k\ + m — 1. 

Likewise, choose k 3 > k 2 , 1 < s < N — n, and consider the new set {r)' k ' : k = 
1, . . . , iV} consisting of the points 

Vkx < ■ ■ • < Vki+s-2 < Vki+s < ■■■ <Vk 2 < Vk 3 

and the corresponding N x iV-matrix Pq = (g(C,j — f]^))j^=i,...,N- This matrix is 
obtained from P by adding the column (g(£j — Uk 3 ))i<j<N as its last (=iV-th) 
column and deleting the column — Vs))^^- Then ?77v-n+i appears in the 

n + 1-st column of Pq, counted from right to left, and 

^N-n — WN-n+l = V^+N-n = Vk 2 -n+l > £n 

by (120 p . Again the Schoenberg- Whitney conditions are violated and therefore 
detPp = 0. We conclude that the vector (g(£j — Vk a ))i<j<N must lie in the linear 
span of the last n columns of Po. 

Step 5: Construction of the left-inverse. Recall that k\ = —m{r + 1) + 1 
and that ?7( r+ i) m = y^+ir+^m-i — Ho- Let ° T denote the (r + l)m-th row vector 
of Po" 1 . By definition of the inverse, we have YljLi c j9(Cj ~ Vk) — 8k,(r+i)m, or 
equivalently, for k\ < k < k 2 , 

N 

(21) ^cjg^j -y k ) = 5 kfi . 

i=i 

Let us now consider the other columns with k < k\ or k > k 2 . Since ki < —m < 
< n < k 2 and every vector — Vk))i<j<N lies either in the span of the first m 
columns of Po (f° r k < ki) or in the span of the last n columns of Po (for k > k 2 ), 
we obtain that 

N 

Yl c ^(0 - y*) = ■ 

3=1 

Therefore, the identity (|2Tj) holds for all k G Z. 

Next we fill the vector c with zeros and define the infinite vector 70 by 

if xj = G X' Q 
otherwise. 

N 

3=1 



7o, 



c .i' 



Then 



^2lo,j9(xj ~Vk) 
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Thus 70 is a row of the left-inverse of P. By construction, 70 has at most N non- 
zero entries. In particular, if Xj < y^+m-i = V-rm, then we have j < ji and thus 
70 j = 0. Similarly, if x, } > y k2 - n +i = y r {n-x) for n > 1, and Xj > y k2 -n+i = V\ for 
n — 1, then we have j > j 2 and 7oj = 0. This gives the support properties of the 
entries 70^ of row k = of the left-inverse T. 

Step 6: The other rows of T. The construction of the k-th row of T is similar. 

We choose the columns between k\ = k — (r + l)m + 1 and k% = k\ + N — 1 of P, 
and, accordingly, we choose suitable rows between the indices ji = min{j : Xj > 
y kl+m -\} and j 2 := max{j : Xj < y k2 -n+i}- Then the column of P containing y k 
has at least m columns to its left and n columns to its right. We then proceed to 
select £j's and define an N x N- matrix P k and verify that det P k > 0. The k-th 
row of T is obtained by padding the appropriate row (with row-index (r + l)m) of 
PfT 1 with zeros. By this construction one obtains a vector r y k = {'Jk > j)jez, for which 

^2lkjg(xj - yi) = 5 kj i 

holds. Furthermore, jk,j = when Xj < yk-rm and when Xj > y k +m- 

Step 7: The remaining cases. The cases where m = or n = are simple 
adaptations of the above steps. For m = 0, n > 2 we choose N = (n — l)(r + 1). 
The indices for the submatrix P k that occurs in the construction of the k-th row 
of T are ki = k, k 2 = k + N — 1, ji = j\{k) = j\ = max{j : Xj < y kl } and 
H = 32(k) = max{j : Xj < yk 2 }- Now we proceed as before. 

We note that Step 4 simplifies a bit. For m = the function g is supported on 
(—00, 0) by [251 P- 339]. Consequently, if k < k\ and j > j%, then Xj — y ko > and 
9(xj ~ Vko) =0. 

Thus the column vectors of P to the left of the submatrix Pq are identically zero 
and no further proof is needed for linear dependence. 

The case n = is similar. It can be reduced to the previous case by a reflection 
x — > —x, which interchanges the role of m and n. 

The special case of m = and n — 1 can be solved by taking N — 2, ki — k, 
fc 2 = k + 1, ji — max{j : Xj < y k }, 32 = min{j : Xj > y k }, and the 2 x 2-matrix 

p = ( 9<KX n ~ Vk) g(%h ~ Vk+i)\ 

k \a( x j2 - Vk) g(xj 2 - Vk+i) J ' 

In this simple case the size of the matrix is independent of the parameter r occurring 
in condition C r . ■ 

3. Gabor Frames with Totally Positive functions 

In this section we prove the main result about Gabor frames. Recall that a totally 
positive function is said to be of finite type, if its two-sided Laplace transform 
factors as = Ce Ss rL=i(l + $v s ) with real numbers 5, 8 V . 

Theorem 9. Assume that g is a totally positive function of finite type and M = 
m + n > 2, where m is the number of positive zeros and n the number of negative 
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zeros of 1 / C(g) . Then J 7 (g) = {(a, (3) G R+ : a(3 < 1}. The Gabor frame Q(g,a, (3) 
possesses a piecewise continuous dual window 7 with compact support in [— — 

a, 7j + a], ^/jere r := [j^J . 

By taking a Fourier transform, we obtain the following corollary. 

Corollary 10. If h(r) = CU™ =1 (1 + 2m8 v r)- 1 for M > 2, then F(h) = {(a, (3) G 
]R^_ : af3 < 1} and Q(h, a, (3) possesses a bandlimited dual window 6 with supp 9 C 

REMARK: We have excluded the case M = 1 for the formulation of the theorem, 
because it corresponds to the one-sided exponential g(t) = e~*XR + (i)- This function 
is discontinuous and therefore is not subject to the Balian-Low principle. Instead, 
we have only F{g) = {(a, (3) G R 2 + : a(3 < 1}. 

Proof of Theorem^ Since g G L 1 (R) is totally positive, it decays exponentially, 
and since M > 2, its Fourier transform g(£) = Ce 27Ti5 ^ Ylt = i( 1 + 27ri ^O~ l decays at 
least like \g(£)\ < C(l+£ 2 )~ 1 . In particular, g is continuous. As a consequence, the 
assumptions of the Balian-Low theorem are satisfied [HE] and J-'(g) C {(a, (3) G 
R 2 + : a(3 < 1}. 

To prove that Q(g,a, (3) is a frame for a (3 < 1, we will construct a family of uni- 
formly bounded left-inverses for the pre-Gramians P(x) of g and then use LemmalH 

Fix x G [0, a] and consider the sequences Xj = x + aj and = kf (3, j,k G Z. We 
first check condition (C r ). By our assumption, we have a < 1/(3 and every interval 
(k/p, (k + l)/f3) contains at least one point x + aj. Every interval (k/(3, (k + r)/f3), 
with r G N, contains at least r + 1 points x + cr/', if r/fi > (r + l)o, i.e., we have 

a(3 

r > -, or equivalently r > 

1 — ap 

Consequently, condition (C r ) is satisfied with r = [ 1 _ 1 a/3 J . By Theorem El each 
pre-Gramian P(x) with entries g(x + aj — k/(3) possesses a left-inverse r(x). 

To apply Lemma [5], we need to show that T(x),x G [0, a], is a uniformly bounded 
set of operators on £ 2 (Z). 

Let Pq(x) be the N x A^-square submatrix constructed in Steps 1 and 2. The 
column indices k\ and k 2 depend only on the type of g, but not on x. The row 
indices j x = ji(x) = min{j : xj > y fel+m _i} = min{j : x + aj > (k x + m - 1)/ (3} 
and j'2 = 32(2) are locally constant in x. Likewise the indices that determine 
which rows j,j± < j < ji of P(x) are contained in Pq(x) are locally constant. 
Consequently, for every x G [0, a] there is a neighborhood U x such that indices 
used for P (y)j k — g(y + £j — k/(3) do not depend on y G U x . Since g is continuous, 
Po(y) is continuous on U x , and since detP (x) > there exists a neighborhood 
Kb C C4, such that detP (y) > detP (^)/2 for y G 14- 

We now cover [0, a] with finitely many neighborhoods V x and obtain that 
det Po(y) > min g det P (x q )/2 = 5 > for all y G [0,a]. Since each entry of 
the inverse matrix Po(y) ~ x can be calculated by Cramer's rule, these entries must 



1 

1 -a(3 
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be bounded by CdetPo(?/)~ 1 < C5~ x with a constant C depending only on \\g\\oo 
and the dimension N of Po{y)- 

By construction (Step 5), the zero-th row 7(0;) = {lo,j{ x ))jez °f the left-inverse 
r(zc) contains at most N < (r + 1)M non-zero entries, namely those of the row 
of Pq(x)^ 1 corresponding to y = 0. We have thus constructed a vector-valued 
function x — > j(x) from [0, a] — > £°°(/L) with the following properties: 

(i) j(x) is piecewise continuous, 

(ii) card (supp j(x)) < (r + 1)M, where according to (II 



— V 771 TTl — T7TI TTt 

supp j(x) = {j : 7oj(a;) ^ 0} C {j : — — < x+aj < — } C {j : — - — 1 < j < — } 

p p ap ap 

(iii) and sup^^ ||7(^)|U = C < 00. 

Consequently, the dual window 7(2) = J2jez d 7o,j( x )X[o,a)(2 ; — oej) correspond- 
ing to T(x) by LemmaHl has compact support on the interval [— — a, y + a], is 
piecewise continuous, and is bounded. In particular, it satisfies the Bessel property 
(see [33] or [H Cor. 6.2.3]). 

We have constructed a dual window for Q(g, a, (3) satisfying the Bessel property. 
By Theorem [3] and Lemma HI Q(g, a, /3) is a Gabor frame. ■ 

3.1. Remarks and Conjectures. In the proof of Theorem [9] we have constructed 
a compactly supported dual window 7 for Q(g, a, (3). This construction is explicit 
and can be realized numerically, because it requires only the inversion of finite 
matrices. To determine the values 7(2; + ja), one has to solve the linear N x N 
system P (x)^f(x) = e for a vector e of the standard basis of M N . 

We observe that the canonical dual window (provided by standard frame theory) 
has better smoothness properties. The regularity theory for the Gabor frame oper- 
ator implies that the canonical dual window 7 decays exponentially and its Fourier 
transform 7°(£) decays like (9(|£|~ M ), where M is the type of g. See [12l[P 



Theorem M raises many new questions. Theorem suggests the natural conjec- 
ture that the frame set of every totally positive continuous function g in L 1 (R) is 
J-'ig) = {(a, (3) G 1R^_ : a(3 < 1}. Our proof is tailored to totally positive functions 
of finite type, most likely the proof of the conjecture will require a different method. 

In a larger context one may speculate about the set Ai of functions such that 
the frame set is exactly F(g) = {(a,f3) G : a(3 < 1}. In other words, when 
is the necessary density condition a(3 < 1 also sufficient for Q(g,a,(3) to be a 
frame? The invariance properties of Gabor frames imply that the class Ai must 
be invariant under time-frequency shifts, dilations, involution, and the Fourier 
transform. Furthermore, if F(g) = {(a,/3) G 1R^_ : a/3 < 1}, then both g and g 
must have infinite support. 

A general method for constructing functions in Ai can be extracted from |22j. 
We write c(£) = ^ fcgZ Cke 2mk ^ for the Fourier series of a sequence (c^) and then 
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define, for a given function go G L 2 (IR), 

C go = {fe L 2 (R) :f=J2 CkdiT k M igo , C) de t\Z), inf(|c(0 d{^)\ > 0} . 

k,iez 

Lemma 11. Let g be a totally positive function of finite type M > 2. If g G C go , 
then the frame set of g is T{g) = {(a, (3) G : a(3 < 1}. 

This trick was used by Janssen and Strohmer in [22]. They showed that the 
hyperbolic secant gi(t) = (e* + e - ') -1 belongs to the set C v for the Gaussian 
window (p(t) = e~ nt and then concluded that J-{g\) = {(a,f3) G : af3 < 1}. 
The general argument is identical. 

Each class C g is completely determined by the zeros of the Zak transform of g. 
Let Zg(x, f ) = J] fcgZ g(t-k)e 2nik ^ be the Zak transform of g. Since Z{T k M ig ) (x, f ) = 
e 27T ^ lx ~ k ^Zg(x,^), every (7 £ C 9o has a Zak transform of the form 

Zg(x,0 = c(-Z)d(x)Zg (x,£) . 

The definition of C go implies that Zg and Zg have the same zero set. If go and h 
are two totally positive functions, then the zero sets of Zg and Zh are different 
in general, therefore Lemma [TT] leads to distinct sets C ga . The zeros of the Zak 
transform seems to be some kind of invariant for the Gabor frame problem, but 
their deeper significance is still mysterious. 

4. Sampling Theorems 

In this section we exploit the connection between Gabor frames and sampling 
theorems and prove new and sharp sampling theorems for shift-invariant spaces. 
Originally shift-invariant spaces were used as a substitute for bandlimited functions 
and were defined as the span of integer translates of a given function g. We refer 
to the survey [3] for the theory of sampling in shift-invariant spaces. We will deal 
with a slightly more general class of spaces that are generated by arbitrary shifts. 

Let Y = (yk)kez be a strictly increasing sequence and consider the quasi shift- 
invariant space 

V Y (g) = {fe L 2 (R) :f = J2 c ^(- - Vk)} ■ 

fcez 

We require that the sequence Y = (y k ) of shift parameters satisfies the conditions 

(22) < q Y = inf(jfc+i - y k ) < sup(y fc+ i - y k ) = Q Y < 00. 

k k 

Such sequences are called quasi- uniform or uniformly discrete. The numbers Qy , qy > 
are the mesh-norm and the separation distance of Y. 

For the norm equivalence ||/|| 2 x ||c|| 2 for / G Vy(g) we need that {g(. — y k : 
y k G Y} is a Riesz basis for Vy{g)- 

Lemma 12. Let g be an arbitrary totally positive function. 

(i) IfY— KL with h > 0, then {g(. — hk), k G Z} is a Riesz basis for Vy{g)- 

(ii) IfY is quasi-uniform, then {g(. — y k ) : k G Z} is a Riesz basis for Vy{g)- 
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Proof, (i) Since g is continuous, does not have any real zeros, and g(£) decays 
at least like C/|£|, every periodization of \g\ 2 is bounded above and below. This 
property is equivalent to the Riesz basis property, e.g. Thm. 7.2.3]. 
Of course, (i) also follows from (ii). 

(ii) For the general case we use Zygmund's inequality [351 Thm. 9.1]: If I is an 
interval of length III > — , then 

to II qy 1 

j |^c fc e 2 ^| 2 ^> A 5 |/|||c|| 2 

for a constant depending only on S > 0. 
If / = Y.k c k9(- ~Vk), then 

11/112 = 11/112= / |5> e " 2W ! 2 \9(T)\ 2 dT 

Jr k 

>inf|£(r)| 2 I \Y,c k e- 2mykT \ 2 di 

k 

> l7|/|A 5 ||c|| 2 . 

Here inf rg / |^(t)| 2 > 0, because g does not have any real zeros by Theorem EJ ■ 

We are interested to derive sampling theorems for generalized shift-invariant 
spaces that are generated by a totally positive function g. Our goal is to construct 
strictly increasing sequences X = (xj) that yield a sampling inequality 

(23) A\\f\\l < £ \f(^)\ 2 < B\\f\\l for all / G V Y (g) 

for some constants A, B > independent of /. Following Landau [24], a set X C R 
that satisfies the norm equivalence (|2"31 is called a set of (stable) sampling for Vy{g)- 
Except for bandlimited functions and B-spline generators only qualitative results 
are known about sets of sampling in shift-invariant spaces. 

We first give an equivalent condition for sets of sampling in Vy(g)- As in LemmaH] 
we obtain the following characterization of sets of sampling in Vy{g)- 

Lemma 13. Let g e L 2 (R), and let Y = (yk)kez C R be a strictly increasing 
sequence. Then a set {xj} C R is a set of sampling for Vy{g), if and only if 
the pre-Gramian P with entries pjk = g(xj — y^) possesses a left-inverse V that is 
bounded on £ 2 (Z). 

The case of uniform sampling in shift-invariant spaces is completely settled by 
the results in Sectional 

Corollary 14. Let g be a totally positive function of finite type M > 2 andY = KL. 
If a < h and x G R is arbitrary, then the set x + aZ* is a set of sampling for Vy{g). 
More precisely, there exist positive constants A,B independent of x, such that 

A\\f\\l < \f(* + a J)\ 2 < B Wf\\l for oil f G V Y (g) . 
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Proof. We proved Theorem [9] by verifying the equivalent condition of Theorem |3l 
namely (J1J stating that 

(24) A||c||| < I ^2c k g(x + ja- hk)\ 2 < B\\c\\\ for all i6l,c6 £ 2 (Z) . 

Since / G Vy(g) is of the form / = J2k c kd(- ~ hk) and ||/|| 2 x ||c|| 2 by Lemma [T2| 
the inequalities are equivalent to the sampling inequality ||/||| x \ f( x + 

aj)\ 2 , and the constants are independent of x by Theorem [3j ■ 

REMARK: The condition a < h is sharp. If a = h, then there exists x G R, such 
that i + aZ is not a set of sampling. This follows immediately from the Balian-Low 
theorem jl]. 

Our methods yield more general sampling theorems. On the one hand, we 
study non-uniform sampling sets, and on the other hand, we may treat quasi shift- 
invariant spaces. The auxiliary characterization of Lemma [13] gives a hint of how 
to proceed. If the sequences (xj) and (y^) satisfy condition (C r ) for some r > 0, 
then by Theorem [5] the pre-Gramian matrix P possesses an algebraic left-inverse. 
To obtain a sampling theorem, we need to impose additional conditions on (xj) 
and (Vk), so that this left-inverse is bounded on £ 2 . 



To verify the boundedness of a matrix, we will apply the following lemma which 
is a direct consequence of Schur's test, see, e.g., [TJJ Lemma 6.2.1]. 

Lemma 15. Assume that A = (djk)j,ke% is a matrix with bounded entries \djk\ < C 
for j, k G Z. Furthermore, assume that there exists a strictly increasing sequence 
(jk)kez °f row indices jk G Z and N G N such that ajk = for \j — jk\ > N. Then 

(25) ||A||^ 2 < (2N-1)C. 

Proof. The conditions give 

jk+N-l 

K 2 := sup |ajjfc| = sup \djk\ < (2N — 1)C. 

For the estimate of the column sums, we define the set 

Nj = {keZ: a jk ^ 0} C {A; G Z : \j - j k \ < N} for j G Z. 
Since (jh) is strictly increasing, Nj has at most (2N — 1) elements, and this gives 

K\ : = sup^^ \a,jk\ = sup \djk\ < (2N — 1)C. 



The assertion now follows from Schur's test. 



In the following we give a sufficient condition for a set X to be a set of sampling 
for V Y (g). 
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Theorem 16. Let g be a totally positive function of finite type M > 2. Let 
Y = {yk)k£Z C R be an increasing quasi-uniform sequence with parameters ay, 
Qy defined in (122]) . Moreover, let (xj)j € % C R be a strictly increasing sequence, 
which satisfies the following conditions: 

' There exist r G N ; e G (0, qv/2), and a quasi-uniform subsequence 
X' C X , such that 

(C r (e)) < (a) every interval (y k + c,yk+i — e) contains at least one point 
Xj G X'/ 

k f&j /or every k £ Z, we have \(y k + e, ?/fc +r — e) D X'\ > r + 1. 
T/ien X is a set of sampling for Vy(g). 

Proof. Step 1. First, we construct a left-inverse of the pre-Gramian P as in the 
proof of Theorem [S] with a small modification. 

We consider only the case m > 1 , n > 1 . For the construction of the row with the 
index k of the left-inverse T, we choose the size N G N for a square submatrix P k 
of P as in Step 1 and the column indices k\ = k — rim + 1) + 1 and ki = hi + N — 1 
as in Step 6. 

To incorporate condition C r (e), we modify the selection of the row indices in 
Step 2 as follows. Assumption C r (e) and our definition of N imply that the interval 
/ = (y^+m-i + e, yk 2 - n +i — e ) contains at least N points Xj G X', where X' is the 
quasi- uniform subset of X in condition C r (e). Define 

ji := min{j : xj G X', Xj > yyt 1+m _i+e}, j 2 := max{j : Xj G X', Xj < y k2 - n+ i-e} , 
then the set 

X k = {xj G X' : ji < j < j 2 } C (y^+m-i + e, 2/fe 2 -n+i - e) 

has at least N elements. We now choose one point Xj G (yi + e, yi + \ — e) D X' 
for each ki + m — 1 < / < k 2 — n and an additional n + m — 1 points X,- G 
(Ukr+m-i + 2/fc 2 -n+i — e) fl X' and obtain a subset 

containing precisely X elements. As before, we set 

and define the quadratic submatrix of P by 

^ = (g(tj - vi))j,i=i,.,N ■ 

The modified construction leads to a stronger version of the Schoenberg- Whitney 
conditions, namely 

(26) £j + e < r/j +m for 1 < j < X — m, + e < for 1 < j < N — n . 

Step 4 remains unchanged, and the column of P with k < ki or k > k 2 are 
linearly dependent on the columns of P k - 

Hence P k is invertible and, by padding the (r + l)m-th row of P- 1 with zeros, 
we obtain the row {jk,j)j€Z with the row index k of the left-inverse Y. 

Step 2. We show that this left inverse V defines a bounded operator on £ 2 (Z). 
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By construction the j-th row (jkj) of T has at most N non-zero entries between 
k\ — k — (r+ l)m + 1 and k\ + N — 1. To apply Lemma [T5| we need to show that 
the entries of Y are uniformly bounded, or equivalently, that the entries of P^ 1 are 
bounded with a bound that does not depend on k. 

We set up a compactness argument similar to the proof of Theorem [HJ 

We begin with the simple observation that 

9(€j - m) = - m) - (vi-vi)), j,i = i,...,N. 

Let S be the A-dimensional simplex 

(27) S = {r = (T 1 ,...,T N )eR N : < n < ■ ■ ■ < r N < (N - 1)Q Y } . 

Although the finite sequences (Cj)i<j<N an d (vi)i<i<n depend on the row index k 
(and we should write ^ and ^ to make the dependence explicit), we always 
have 

0<£i-Vi<£n-Vi<Vn-Vi<( n - 1 )Qy ■ 

Consequently, 

(£1 - T]i, . . . ,£jv - 771) e S, and (0, T] 2 - 771, . . . , r) N - 771) G S . 

Let q := min{gx';QY} > be the minimum of the separation distances of the 
quasi-uniform sets X' and Y and let 

S q = {r = (ri, . . . , Tjv) e S : r j+1 - t 3 > q for 1 < j < N - 1} . 
Then S g is compact and 

(^1 - 771, . . . , £iv - 771) G 5 9 and (0, 772 - . . . , r] N - G 5, . 
Finally, we define the compact set 

K = {(t,6) E S q X S q : tj + e < 6 j+m for 1 < j < N - m, 

0j + £< T j+n for 1 < j < N - n } . 

The assumption C r (e) implies that 

((6 - »7i, • • • > ^ - 771), (0, 772 - ?7i> • • • , m - vi)) e ^ • 

Clearly, the Schoenberg- Whitney conditions are satisfied for every point (r, 9) G if 
and therefore every iV x A^-matrix (^(r,- — 6ij) has positive determinant. Since the 
determinant depends continuously on (r, 6) and if is compact, we conclude that 

This construction implies that det > 5 > for every k. As in the proof of 
Theorem [8] we use Cramer's rule and conclude that all entries of P^ 1 are bounded 
by (iV-l)!^ 1 !!^- 1 . 

The assumptions of the modified Schur test are satisfied, and Lemma [15] yields 
that the matrix T is bounded as an operator on £ 2 (Z). Finally, Lemma [T3l implies 
that A is a set of sampling for Vy(g)- I 
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Corollary 17. Assume that g is totally positive of finite order M > 2 and Y = KL. 
Let a = s\ij)j eZ (Xj+i—Xj) be the maximum gap between consecutive sampling points. 
If a < h, then (xj) is a set of sampling for Vy(g)- 

Proof. The assumption of Theorem [TB] is verified with e = h — a. ■ 

REMARKS: 1. So far the conclusion of Corollary [TTl was known only for 5-splines 
as generators [2]. For other generators only qualitative results were known [TJ or 
weak estimates far from the correct sampling density [Ijl3|l3l]. 

2. If Y = hZ and X satisfies condition C r (e), then the largest possible gap of 
consecutive points in X is 2h — 2e. Of course, large gaps have to be compensated 
by a higher density of neighboring points so that the condition \ {y k + e, y^ +r — e) n 
X'| > r + 1 is still satisfied. By refining the compactness argument in the proof 
of Theorem [TH] it is possible to derive even weaker conditions for X to be a set of 
sampling for a shift-invariant space Vy(g) with totally positive generator g. 
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